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SUMMARY 

The  polar  orbit  of  Heos  2  second-stage  rocket,  1972-05B,  has  been  determined  on  each 
of  the  final  16  days  before  its  decay  in  September  1978,  using  the  RAE  orbit  refinement 
program,  PROP  6,  with  abovt  1360  observations.  An  accuracy  of  30-70  m,  both  radial  and 
across  track,  was  achieved. 

Eleven  values  of  density  scale  height  have  been  determined  from  the  decrease  in 
perigee  height,  with  a  2  per  cent  error;  seven  of  these  values  are  within  6  per  cent  of 
the  CIRA  1972  reference-atmosphere  values,  the  rms  value  being  4  per  cent  higher  than 
CIRA. 

The  rotation  rate  of  the  upper  atmosphere,  A  ,  was  determined  from  the  decrease  in 
orbital  inclination  as  A  •  1.40  ±  0.05  rev/day;  ie  a  strong  west-to-east  zonal  wind  of 
160  ±  20  m/s,  at  a  mean  height  of  about  240  km.  The  local  time  was  01-02  h;  solar 
activity  was  high;  and  the  latitude  of  perigee  moved  steadily  from  10°N  to  67°S. 
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1  INTRODUCTION 

The  ESRO  highly-eccentric-orbit  satellite  Heos  2  (1972-05A)  was  launched  into  an 
orbit  of  eccentricity  0.95  on  1972  January  31  by  NASA.  The  second-stage  Ablestar  rocket, 
1972-05B,  entered  a  much  lower  eccentricity  polar  orbit  (e  =  0.06)  with  a  life  of  about 
6j  years  and  numerous  orbits  are  being  determined  at  the  University  of  Aston;  here  the 
orbit  is  analysed  only  in  the  final  16  days  of  its  orbital  life,  1978  September  12-27. 

The  North  American  Air  Defense  Command  (NORAD)  kindly  supplied  about  2500  observations 
for  this  16-day  period. 

The  RAE  orbit  refinement  program*,  in  the  PROP  6  version  (which  has  a  limit  of 
100  observations  per  orbit),  has  been  used  to  determine  16  daily  orbits  with  about  1500  of 
the  NORAD  observations,  and  an  additional  70  from  the  US  Navy  Navspasur  system. 

The  upper-atmosphere  rotation  rate  has  been  determined  from  the  decrease  in  orbital 
inclination  to  give  the  zonal  wind  speed  near  240  km  height,  and  the  density  scale  height 
has  been  evaluated  over  the  height  band  200-270  km  from  the  decrease  in  perigee  height. 

2  THE  ORBIT 

About  1360  observations  were  used,  with  a  further  200  or  so  rejected,  in  the  final 
determination  of  the  16  daily  orbits  for  1978  September  12-27,  the  last  orbit  being  at  an 
epoch  13  hours  before  decay.  The  orbital  elements  derived  are  listed  in  Table  1,  with 
standard  deviations. 

The  sd  in  inclination,  i  ,  varies  between  0.0003  and  0.0006°,  and  in  eccentricity, 
e  ,  from  4  x  10^  to  7  x  10^.  These  are  equivalent  to  between  30  and  70  m  in  position 
for  both  i  and  e  .  This  is  an  excellent  accuracy  for  an  orbit  so  near  decay.  The  sd 
in  right  ascension  of  the  node,  0  ,  is  between  0.0002  and  0.0008°;  for  argument  of 
perigee,  lo  ,  or  mean  anomaly  at  epoch,  ,  the  sd  varies  from  0.03  to  0.19°.  Figs  1 
to  5  show  the  variations  of  i,  a,  e,  0  and  u  .  Although  the  overall  drop  in  inclina¬ 
tion  over  the  16  days  was  only  0.02°  (see  Fig  1),  it  was  sufficient  to  determine 
accurately  the  upper-atmosphere  wind  (section  5). 

3  VARIATION  IN  PERIGEE  HEIGHT 

Perigee  height  is  required  for  determining  density  scale  height  and  the  heights 
applicable  to  density  scale  height  and  atmospheric  rotation.  Values  of  semi  zaajor  axis, 
a  ,  and  eccentricity,  e  ,  from  Table  1  are  used  to  calculate  perigee  height  over  a 
spherical  Earth,  h  ,  where 


h«a(l-e)-R,  (1) 

P 

where  R  ,  the  Earth's  mean  equatorial  radius,  is  taken  to  be  6378.14  km.  Next,  the 

perturbation  Ae  ,  in  e  ,  due  to  zonal  harmonic  and  lunisolar  effects,  is  obtained 

2 

using  the  PROD  computer  program  and  cleared  from  h  to  give  Q  where 


Q 


h  ♦  aAe  , 
P 


(2) 


lA 

_ 

«n 

CM 

NO 

CM 

o 

o 

CO 

O 

00 

ON 

O' 

00 

ON 

00 

fs. 

O' 

ON 

nO 

c^ 

ON 

ON 

ON 

ON 

ON 

NO 

o 

ON 

On 

ON 

O' 

00 

O' 

ON 

ON 

O' 

o 

d 

o 

d 

o 

O 

— 

o 

O 

o 

d 

d 

o 

d 

d 

d 

ON 

CM 

CM 

00 

_ 

CM 

ON 

O' 

On 

NO 

o 

h) 

m 

CM 

CM 

cn 

CM 

CM 

CM 

CO 

CO 

CO 

CO 

4 

4 

o 

o 

o 

d 

o 

d 

d 

o 

d 

d 

o 

d 

d 

d 

O 

d 

« 

CO  — 

in  CM 

£ 

CM 

NO 

rs  CM 

O  CM 

m  CM 

CO  CO 

CM 

CM 

r^ 

r*. 

OO  CM 

£ 

O 

d 

d 

NO 

OO  vO 

00  't 

—  m 

o  m 

00  — 

nO  *0 

CO  CO 

no  m 

m  m 

■4  m 

OO  00 

On  m 

CO  rs 

CN 

in 

CM 

ON 

rs. 

CM 

o  — 

m 

•4“ 

rs. 

00 

CO 

m 

CM 

NO 

in  CO 

o 

o 

CM 

•o 

00 

o 

fs. 

CM 

CM 

£ 

• 

• 

• 

• 

4 

■ 

, 

<N 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CO 

fO 

m 

NO 

rv 

CM 

d 

CO 

CM  — 

— 

On  — 

00  — 

rs 

_ 

NO  m 

m  CM 

_ 

NO  -s 

in  CO 

CM  — 

CM  CM 

fO  cn 

fs.  CO 

m 

00 

CM 

o 

mm 

o 

NO 

«4 

NO 

o 

rs 

m 

00 

m 

NO 

ON 

CM 

m 

CM 

CM 

•4 

O' 

o 

s4 

s4 

lA 

o 

CM 

m 

£ 

m 

rv 

... 

NO 

o 

m 

d 

m 

d 

rs 

-4 

fO 

sj 

00 

m 

m 

NO 

NO 

rs 

Cv 

00 

00 

O' 

ON 

O 

» 

CM 

CO 

m 

ON 

r>* 

r*«. 

r*^ 

rs 

^•4 

r-s 

rs 

fs, 

rs 

CO 

00 

00 

CD 

OO 

00 

m 

in 

in 

m 

m 

m 

in 

m 

in 

in 

m 

m 

m 

m 

m 

m 

sO 

m  cn 

—  cn 

ON 

CM 

cn 

in  *0 

<4  CO 

fs.  CO 

fs.  m 

fo  rs 

CM  <0 

O  ON 

On  — 

4  >0 

OO  ON 

NO 

00 

nO 

m 

ON 

cn 

CM 

NO 

CO 

o 

fO 

in 

fs. 

On  — 

4 

O'  — 

£ 

d 

m 

m 

ON 

rs 

ON 

CM 

NO 

fs. 

r. 

fs. 

On 

00 

ON 

00 

00 

00 

ON 

CM 

m 

00 

m 

O 

fs. 

nO 

00 

CM 

CM 

fO 

CM 

sO  m 

CM  n 

CM  cn 

m 

NO 

cn 

nO  *0 

O  Mf 

rs  fo 

CO  -4 

CM  m 

CO  rs 

CO  NO 

CO  O' 

in  — 

4  nO 

ON  O' 

nO 

NO 

CM 

On 

o 

m 

m 

NO 

m 

m 

nO 

00 

—  "• 

m 

'O  — 

3 

ON 

in 

r*. 

CM 

00 

00 

CO 

CO 

CO 

fs. 

8 

d 

CM 

in 

lA 

CO 

CO 

CM 

<M 

mm 

o 

O' 

cn 

cn 

cn 

cn 

CO 

CO 

CO 

CO 

fO 

CO 

CO 

fO 

CM 

—  m 

nO 

ST  o.t 

O 

ON 

Mr 

cn  *0 

sj  CO 

CM  CM 

s»  *4 

NO  NO 

fs.  m 

fO  s4 

O  CO 

nO  00 

m  4 

00 

hN 

05 

00 

r*4 

00 

rs. 

rs 

NO 

*4 

rs 

NO 

P 

O' 

r*«. 

cn 

O' 

in 

fs» 

CO 

O' 

m 

fs. 

fO 

4 

O 

m 

a 

NO 

NO 

m 

m 

m 

Mf 

S» 

CO 

CO 

CO 

CM 

CM 

o 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

OO 

00 

00 

—  in 

NO  'O 

o  in 

00  in 

CN.  4.^ 

00  in 

in  CO 

CM  *4 

ON  >4 

O'  in 

vo  in 

—  *4 

in  *4 

rs  \0 

in  in 

r* 

NO 

NO 

n 

CO 

CM 

CO 

OO 

fs. 

nO 

4 

CM 

fs, 

cn 

n 

cn 

rN 

cn 

CO 

CO 

CO 

CO 

CM 

CM 

CM 

CM 

CM 

rv 

r> 

fs. 

rv. 

rs 

rs 

rs 

rs 

r^ 

fs 

rs 

r. 

rs 

fs. 

ON 

ON 

O' 

ON 

O' 

O' 

O' 

O' 

O' 

ON 

O' 

O' 

O' 

ON 

ON 

O' 

■ 

00 

00 

00 

00 

00 

00 

00 

00 

00 

OO 

00 

CO 

00 

00 

CO 

00 

rn  vo 

O  MT 

—  m 

O  m 

ON 

m 

CM 

00  rs 

00  «4 

o  m 

\0  in 

rs  so 

m  m 

fs.  m 

4  4 

4  rs 

m  no 

cn 

O 

ON 

00 

rs. 

NO 

rs 

ON 

00 

mm 

O' 

OO 

fO 

mm 

O 

p. 

<N 

O' 

in 

CM 

O' 

NO 

m 

O 

rs 

m 

CO 

in 

mm 

«0 

00 

fs. 

NO 

NO 

nO 

m 

in 

in 

m 

S^^ 

s4 

*4 

CO 

fO 

fO 

CM 

V 

o 

g 

o 

O 

p 

O 

o 

o 

8 

o 

O 

O 

O 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

O 

O 

o 

o 

o 

o 

d 

d 

d 

o 

o 

d 

d 

o 

d 

d 

d 

d 

d 

d 

d 

d 

—  o 

CM  O' 

cn  O 

CM 

00 

\0  rs. 

O  CO 

00  00 

NO  — 

—  CO 

m  O' 

nO  s4 

tn  no 

CM  O 

00  O 

\0  fs. 

CM  — 

NO  — 

ON 

ON  — 

NO 

O'  n 

sj  — 

rs 

00  — 

m  CM 

CM 

nO  — 

CM  CM 

CO  CM 

O  4 

in  \o 

rv 

CM 

00 

CM 

rN 

o 

CO 

00 

m 

CM 

rs 

o 

CM 

rs 

m 

o 

NO 

in 

rs 

rs 

-4 

OO 

CM 

*4 

O' 

00 

«£ 

00 

On 

NO 

CM 

ON 

NO 

CM 

00 

•4 

o 

m 

o 

CO 

s4 

4 

ON 

m 

m 

in 

CO 

CO 

CO 

CM 

CM 

O 

O' 

fs. 

m 

nO 

NO 

NO 

NO 

NO 

NO 

NO 

NO 

yO 

NO 

nO 

NO 

nO 

in 

m 

m 

nO 

NO 

nO 

NO 

NO 

NO 

nO 

NO 

NO 

NO 

nO 

NO 

NO 

nO 

NO 

o 

o 

O 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 

p 

m 

nO 

r«« 

00 

O' 

o 

CM 

CO 

s4 

in 

>0 

fs. 

00 

NO 

nO 

NO 

NO 

NO 

NO 

NO 

rs 

rs 

fs. 

rs 

fs. 

fs. 

fs. 

rs 

J* 

r*» 

rN 

«.» 

■ 

a 

2  w 

CM 

O 

•o 

m 

NO 

rs. 

00 

O' 

O 

CM 

fO 

s4 

m 

nO 

rs 

2s 

“ 

“ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

ff' 

w 

■■ 

jm 

.«4 

mm 

CM 

m 

00 

ON 

O 

mm 

CM 

4 

in 

NO 

A 

wm 

** 

o 

■1 

J 


which  should  be  free  of  gravitational  perturbations.  These  values  of  Q  are  plotted  in 
Fig  6,  and  its  decrease  iQ  is  used  in  the  determination  of  scale  height  (section  4). 


Perigee  height,  ,  over  an  oblate  Earth  is  obtained  by  allowing  for  (a)  the 
variation  in  local  Earth  radius  at  perigee  latitude  and  (b)  the  dynamical  effects  of  the 

3 

Earth's  oblateness  ,  to  give 

y  =  h  +  21.38  sin^o)  +  (1  .67  -  3.32  sin^u)  , 

P  ’ 


for  1972-05B  at  inclination  90  .  Hence 


y  =  h  +  18.06  sin^o)  +  1.67  km  (3) 

P  P 

and  is  also  plotted  in  Fig  6. 

4  DENSITY  SCALE  HEIGHT 

The  decrease  in  Q  ,  shown  in  Fig  6  to  be  very  well  defined,  depends  on  the  density 
scale  height  H  ,  and  analysis  of  the  variation  in  Q  should  therefore  yield  reliable 
values  of  H  . 

Lox^eccentricity  theory  is  needed  because  z  (=  ae/H)  is  less  than  3,  and,  for  an 

.  .  4 

orbit  in  an  oblate  atmosphere,  we  have 


da 

dx 


ie(4  -  3y^  -  y^y^)  -  Jc  cos  2«(yQ  -  2y^  +  y^y^) 


(4) 


where  x  =  ae  ,  y^  =  >  ij.  being  the  Bessel  function  of  the  first  kind  and 

imaginary  argument,  of  order  r  and  argument  z  ;  and  c  =  c'a(l  -  e)  sin^i/2H  ,  with 
e'  ,  the  ellipticity  of  the  atmosphere,  taken  equal  to  the  Earth's  ellipticity,  0.00335. 

In  previous  orbital  analyses  using  low-eccentricity  theory,  the  perigee  height  has 

been  below  200  km  and  equation  (4)  has  been  used  as  it  stands.  Here,  however,  values  of 

H  are  being  determined  at  heights  up  to  about  260  km,  where  the  day-to-night  variation 

in  air  density  is  large  enough  to  affect  equation  (4)  appreciably.  In  an  atmosphere  with 

maximum  daytime  density  p  (at  14  h  local  time)  and  minimum  night-time  density  p  . 

max  ®  j  ■'  '^min 

(at  02  h  local  time),  the  equation  for  da/dx  for  a  low-eccentricity  orbit  is 


da 

dx 


_ Yq  ^  F  cos 

I  +  \F  cos  *  y2^ 


{1  +  0(e)) 


(5) 


where  is  the  angular  distance  of  perigee  from  the  centre  of  the  diurnal  bulge, 

F"  (f-  l)/(f  +  1)  and  f  ”  p  /p  .  .  In  the  conditions  experienced  by  1972-05B 

(exospheric  temperature  at  night  =  900  K) ,  the  COSPAR  International  Reference  Atmoaphere 
1972  (Ref  6)  gives  f  »  1.3  at  200  km  increasing  to  f  -  1.7  at  260  km  height,  so  that 
0.15  <  F  <  0.25.  Also  z  ~  I  ,  so  that  j^y^  *  ~  *  ■  Assuming  that  equation  (5)  can 

be  expanded  in  powers  of  F  ,  we  find 


6 


da 

dx 


(6) 


If  we  assume  that  the  effects  of  c  and  F  are  simply  additive,  the  equation  for  da/dx 
in  an  oblate  atmosphere  with  day-to-night  variation  in  density  may,  from  (4)  and  (6),  be 
wri tten 


^  =  y^  Me(4  -  3yJ  -  y^yj  -  }c  cos  2u>(y^  -  2y,  +  y^y^) 

-  jF  cos  <'p(yQ  ♦  yQy2  "  2)  0(L",F^,e‘,etc)  .  (7) 

For  I972-05B  the  local  time  at  perigee  is  near  02  h,  diametrically  opposite  the  centre  of 
the  diurnal  bulge  when  ei  =  0  ,  so  that,  by  a  convenient  chance,  =  |  1 80°  -  u(  .  Thus, 

for  I972-05B  in  its  final  days,  equation  (7)  becomes 


g  =  yo  +  Je(4  -  3yJ  -  y^yj  -  jc  cos  2u(y^  -  2y ^  *  y^y^) 

*  jF  cos  lij^yp  +  y^y^,  -  2^  =  B  ,  say  .  (8) 

This  theoretical  value  of  da/dx  must  be  matched  with  the  observational  value 
which,  following  Ref  7,  is  given  by 


da 

dx 


.  3AQ  *  ”.b\ 


a  ,  say 


(9) 


where  AQ  is  the  change  in  Q  over  a  time  interval  At  ,  selected  to  ensure  that 

'.'.q!  ^  3  km  ,  and  M,.  and  M,_  are  the  initial  and  final  values  of  the  mean  motion  M. 

1 A  I B  ) 

during  the  time  interval  At  .  In  (8)  and  (9)  the  values  of  a,  e,  cos  2ui  and  cos  a) 
are  taken  as  the  mean  during  this  time  interval. 

The  density  scale  height  H  can  now  be  determined  by  guessing  two  values  of  H  , 
say  H|  and  ,  which  give  Sj  and  62  from  equation  (8),  and  then  by  interpolation 

using 


H 


H|(B2  -  a)  +  H^ia  -  6,) 


(10) 


These  observational  values  of  H  are  plotted  in  Fig  7  as  circles  with  vertical  rms  error 
bars  (calculated  to  be  near  I  km  for  all  values  of  H  )  and  horizontal  bars  representing 
the  time  intervals  over  which  the  values  of  H  are  calculated. 

Values  of  H  have  also  been  calculated  using  the  CIHA  1972  reference  atmosphere^, 
for  the  appropriate  heights  and  exospheric  temperatures,  T^  ,  and  are  plotted  in  Fig  7 
as  crosses.  The  values  of  T^  are  calculated  using  appropriate  values  of  solar  10.7  cm 


radiation  energy  (plotted  in  Fig  7)  and  geomagnetic  index,  with  allowance  for  semi-annual 

g 

effects  .  Care  must  be  taken  here  to  ensure  that  local  time  used  to  calculate  T  is  not 

00 

that  at  perigee  but  averaged  over  a  wide  arc  of  the  orbit,  as  discussed  in  Ref  7.  Using 
local  time  at  perigee  would  result  in  errors  of  up  to  100  K  in  .  Also  plotted  in 
Fig  7  is  the  height  y  (=  yp'*'Y^p)  which  H  applies.  Values  of  y  ,  obtained  from 
Ref  9,  vary  between  0.4  and  0.8. 

Of  the  eleven  observational  values  of  H  calculated,  with  a  2  per  cent  error,  seven 

are  within  6  per  cent  of  the  CIRA  1972  values,  the  rms  value  being  4  per  cent  higher  than 

CIRA.  This  agrees  with  a  previously'  determined  value  of  7  per  cent  higher,  also  for 
high  solar  activity.  The  four  remaining  values  are  10-15  per  cent  higher  than  CIRA  and 
Fig  7  shows  this  occurs  during  a  period  of  high  solar  radiation  energy.  A  medium  magnetic 
storm  occurred  during  the  final  3  days  in  orbit  (1978  September  25-27),  where  Ap 
increased  from  36  to  51,  but  this  does  not  seem  to  have  made  any  significant  difference. 

The  observational  values  of  density  scale  height,  H  ,  plotted  in  Fig  7  are 

replotted  in  Fig  8  against  height  y  ,  as  crosses  against  a  background  of  density  scale 
height  curves,  for  various  T^  ,  obtained  from  CIRA  1972.  Adjacent  to  each  cross  is  the 
value  of  T^  ,  in  brackets,  which  was  used  to  calculate  the  CIRA  value  of  H  . 

Both  atmospheric  oblateness  e'  and  the  day-to-night  variation  in  density, 
measured  by  F  ,  strongly  influence  the  values  of  H  obtained.  If  c  is  taken  as  zero  in 
equation  (8),  ie  for  zero  e'  ,  the  value  of  H  obtained  is  reduced  by  20  per  cent  in  the 
first  few  days,  is  unchanged  on  September  23-24  (when  cos  2a)  =  0  ),  and  is  increased  by 
20  per  cent  on  the  day  of  decay.  If  F  is  taken  as  zero  in  equation  (8),  H  increases 
by  30  per  cent  except  for  the  final  4  days  (to  decay)  when  the  increase  is  about  20  per 
cent.  There  is  no  reason  to  suspect  that  the  values  of  e'  or  F  used  are  appreciably 
in  error;  but  if  they  were  in  error  by  10  per  cent,  say,  errors  of  up  to  2J  per  cent  in 
H  would  occur. 

5  UPPER-ATMOSPHERE  ZONAL  WIND 
5. 1  Results 

The  16  daily  values  of  inclination  determined  are  listed  in  Table  I  and  plotted  in 

Fig  1 .  The  values  were  then  cleared  of  lunisolar  and  geopotential  perturbations  using  the 
2 

PROD  program  ,  with  numerical  integration  at  1-day  intervals.  The  modified  values  are 
plotted  in  Fig  9,  with  standard  deviations.  The  overall  decrease  in  inclination  is  about 
0.023°,  about  50  times  the  mean  sd  of  the  values. 

Since  the  orbital  inclination  is  near  90°,  the  meridional  wind  has  no  effect  and  so 
the  variation  in  inclination  is  due  to  the  zonal  wind. 

The  theoretical  change  in  inclination  has  been  calculated  for  several  values  of 
atmospheric  rotation  rate  (expressed  as  A  times  the  Earth's  rotation  rate)  using  the  SAE 
computer  program  ROTATM  (for  an  oblate  atmosphere),  with  daily  integration  steps.  The 
best  fit  to  the  modified  inclination  values  is  given  by  A  •  1.40  1  0.05  rev /day  and  is 
shown  in  Fig  9,  where  most  of  the  points  are  within  2  sd  of  the  curve.  This  corresponds 
to  a  zonal  west-to-east  wind  of  160  ^  20  m/s.  Since  zonal  wind  speed  depends  on  local 
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time,  this  is  plotted  along  the  top  of  Fig  9;  it  shows  little  variation,  always  being 
between  01  and  02  h.  The  orbits  are  determined  for  1978  September,  ie  equinox  and  high 
solar  activity  conditions;  and  the  latitude  of  perigee  steadily  changes  from  I0°N  to  67°S. 
So,  although  the  effect  of  aerodynamic  forces  is  spread  over  a  wide  arc  around  perigee, 
there  is  a  bias  towards  equatorial  rather  than  polar  regions. 

The  wind  speed  is  also  an  average  over  a  range  of  heights,  but  the  mean  height  y^ 
is  given  by 

y  =  y  +  XH 
A  •'p  p 

where  =  0.5  ,  from  Ref  II.  The  mean  perigee  height  y  is  estimated  to  be  near 

220  km  from  Fig  6,  while  the  mean  perigee  scale  height  H  is  estimated  to  be  near 

40  km.  Hence  y,  =  240  km  . 

■'A 

5 . 2  Discussion 

1  2 

The  analysis  of  numerous  orbits  by  King-Hele  and  Walker  has  indicated  that  zonal 
winds  at  heights  of  200-400  km  are  usually  from  west  to  east  at  local  times  from  18  to  24  h, 
and  east  to  west  from  04  to  12h local  time.  So  the  west-to-east  wind  of  160  m/s  at  240  km 
height  at  01-02  h  local  time  obtained  here  is  unexpectedly  strong.  However,  it  is  not  so 
anomalous  as  appears  at  first  sight,  for  three  reasons. 

(1)  The  zonal  winds  are  luora  strongly  from  west  to  east  near  the  equator  than  at  higher 
latitudes  (see  Fig  14  of  Kef  12),  and  the  perigee  of  I972-05B  travelled  from  latitude 
I0°N  to  latitude  near  60°S,  so  that  the  effective  average  latitude  was  probably  about  20°. 

(2)  The  variation  of  wind  speed  with  local  time  depends  on  latitude,  and  the  maximum 
west-to-east  wind  probably  occurs  at  a  later  local  time  near  the  equator.  For  example, 
the  airglow  measurements  of  Sipler  and  Biondi'^,  for  20  August  1977  at  latitude  9°N, 
indicate  a  west-to-east  wind  exceeding  100  ra/s  between  21  and  02  h,  with  a  value  of 

125  m/ s  at  01  h. 

(3)  The  variation  with  local  time  may  also  depend  on  solar  activity,  with  the  maximum 

west-to-east  wind  being  somewhat  later  when  solar  activity  is  high.  For  example,  in  1970, 

1  4 

when  solar  activity  was  high,  Roble,  Salah  and  Emery  found  that,  at  240  km  height  at 
latitude  42°N,  zonal  winds  were  west  to  east  with  velocity  exceeding  100  m/s  from  18  to 
01  h  on  23-24  February  and  from  19  to  01  h  on  lri-19  May.  But  in  1974,  when  solar  activity 
was  low,  Hernandez  and  Roble found  west-to-cast  winds  exceeding  100  m/s  from  19  to  23  h 
at  240  km  on  I  5  August. 

The  high  value  of  A  obtained  here,  1.40  t  0.05  rev/day,  is  paralleled  by  results 
from  three  other  satellites,  I969-20B,  I969-94B  and  I970-43B,  which  all  give 
A  =  1.3  *  0.1  for  heights  near  240  km  at  time®  of  high  solar  activity,  for  local  times 
of  18-03  h  (see  Ref  12). 

The  zonal  wind  speed  in  the  thermosphere  probably  varies  with  height,  season, 
latitude,  solar  activity,  local  time  and  day  of  the  year,  as  well  as  being  correlated 
with  geomagnetic  disturbances  and  events  in  the  lower  atmosphere.  The  value  from  I972-05B, 
being  so  well  defined,  should  help  in  unravelling  this  complex  pattern. 
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CONCLUSIONS 


The  near-polar  orbit  of  Heos  2  second-stage  rocket,  I972-05B,  has  been  determined 
daily  for  the  last  16  days  in  orbit,  including  an  orbit  for  the  day  of  decay 
1978  September  27.  About  1360  observations  were  used,  most  of  which  were  kindly  supplied 
by  the  North  American  Air  Defense  Command,  NORAD. 

For  the  16  orbits  obtained,  given  in  Table  1,  the  sd  in  inclination  varies  between 
0.0003  and  0.0006°,  while  the  sd  in  eccentricity  varies  from  4  x  10  to  7  x  10  these 
variations  are  both  equivalent  to  30-70  m  in  position. 

Values  of  density  scale  height  have  been  determined  from  the  decrease  in  perigee 
height  and  compared  with  the  CIRA  1972  reference  atmosphere;  the  values  obtained  are 
probably  accurate  to  about  2  per  cent  and  are  about  4  per  cent  above  CIRA,  on  average. 

A  value  for  the  atmospheric  rotation  rate  A  was  determined  from  the  decrease  in 
orbital  inclination  and  found  to  be  A  =  1.40  i  0.05  rev/day  ,  equivalent  to  a  strong 
west-to-east  zonal  wind  of  160  ±  20  m/s  at  a  mean  height  of  about  240  km.  The  conditions 
prevailing  were:  local  time  01-02  h;  equinox;  high  solar  activity;  and  latitude  of 
perigee  between  10°N  and  67°S. 
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Fig  8  Values  of  density  scale  height,  with  CIRA  1972  curvet 
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